Extended states in disordered systems: role of off-diagonal correlations 
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We study one-dimensional systems with random diagonal disorder but off-diagonal short-range 
correlations imposed by structural constraints. We find that these correlations generate effective 
conduction channels for finite systems. At a certain golden correlation condition for the hopping 
amplitudes, we find an extended state for an infinite system. Our model has important implica- 
tions to charge transport in DNA molecules, and a possible set of experiments in semiconductor 
superlattices is proposed to verify our most interesting theoretical predictions. 
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Electronic states in disordered system have been an 
active research topic for many years. Ever since the pi- 
oneering work of Andersonfi it has been generally be- 
lieved that disorder in low dimensional systems leads to 
unequivocal localization of electrons. However, the situ- 
ation changes if additional structure or correlations are 
imposed on the statistical properties of the randomness. 
It was found, for example, that a few special extended 
states in a ID "random dimer" model exist due to sym- 
metries of the resonant scattering in the structured Fur- 
thermore, the existence of a mobility edge separating ex- 
tended and localized states was confirmed for ID random 
systems with weak long-range correlated disorder^ 

Most studies have concentrated on diagonal disorder, 
where the local energies in a tight-binding description are 
assigned randomly, although some studies have explored 
the role of off-diagonal disorder, where the intersite hop- 
ping constants are chosen from a random distribution. 
The role of correlated diagonal and off-diagonal disorder 
has received attention only recently, both theoretically 4,5 
and experimentally 6 Moreover, in many systems, local 
correlations appear naturally due to the built-in chem- 
ical structure. In this article, we investigate the effects 
of structural constraints on the correlated diagonal and 
off-diagonal disorder, and their impact on charge trans- 
port. We find that the local correlations generate ex- 
tended states, which therefore enhance electronic trans- 
port even in the macroscopic limit. 

Our studies have been motivated in part by questions 
on the nature of charge transport in DNA, a subject 
which has arisen much interest recently, due to its funda- 
mental roles in biological processes and in possible novel 
device designs. 7 A DNA molecular system can be viewed 
as a f D chain composed of base pairs AT and CG in a 
typically random order. The on-site energies for pairs 
of bases AT and CG are different, corresponding to the 
different ionization potentials A The role of onsite energy 
correlations in different DNA sequences was discussed re- 
cently in the literature^ where a constant hopping am- 
plitude was considered. However, the hopping ampli- 
tude (via 7r-orbital overlap) in DNA molecules depends 
on whether the electron (or hole) hops between AT/ AT, 



AT/CG, or CG/CG base pairs^ The short-range corre- 
lation of the hopping amplitudes due to the built in chem- 
ical structure is shown to affect the transport properties 
and effectively open conduction channels for electrons in 
DNA molecules - even in those with fully random se- 
quences, such as A-phage DNA. The transport properties 
are shown to be actually determined by a subtle compe- 
tition between the disorder in base pair arrangement (on- 
site disorder) and hopping ("off-diagonal") correlations. 

The minimal model to study random systems with di- 
agonal and off-diagonal disorder is an effective ID tight- 
binding model described by the Hamiltonian 

H = T,j[ £ 3 c t C 3 + *w+i(4 c i+ 1 + 4+i c j)L (!) 

where the onsite energies are chosen from the bivalued 
distribution Sj = ea and Eb- Correspondingly, the hop- 
ping constants are given by tjj+i = tAA (or tss), if 
Ej = £j+i = £a (or Ebb)', while tjj+i = £ab, otherwise. 
This model is perhaps the simplest generalization of the 
Anderson model, which is the limit for tAA = tBB = ^ab ■ 
Notice that in DNA, the A and B labels refer to the two 
kinds of base pairs, AT and CG, while the model could 
be easily adapted to describe electronic states in other 
complex molecules (polymers) and/or semiconductor su- 
perlattices, as we will discuss belowi 1 ^ 

When the concentration of one type of site is small, 
say B, the probability for two nearby sites to have the 
same onsite energy eb is smaller. In this case, the system 
tends to the "repulsive binary alloy" model, in which one 
extended state exists^ A simple calculation yields the 
transmission coefficient for a system with one impurity 
with onsite energy eb, 

Tl{E) ~ {2ii B sinkf + Nr ( ) 

where N\ — Wt A A + 2(t AB — t 2 AA ) cos k, E — Haa cos k, 
and W = £b — £a- One can see that for the state with en- 
ergy E — Wt AA / (t AA — ^ab)' the transmission coefficient 
is unity. The states near this energy have large transmis- 
sion coefficient and long localization length, even in sys- 
tems with more B impurities. In fact, these states have 
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an important contribution to transport. Figure 1 shows 
the transmission coefficient for various concentration of B 
impurities with energy e B . The transmission is obtained 
by a transfer matrix calculation for 1000 sites, and aver- 
aged over 300 different configurations. For the purpose of 
comparison, we also show the transmission coefficient for 
the Anderson model (where t A A —^bb = ^ab = t) with 
the same degree of onsite disorder. We see that in all 
cases, the local correlation built-in through the t values 
leads to much larger transmission coefficient, compared 
with those in the Anderson model. When the concentra- 
tion is small, there is in fact a regime of high transmission 
(~ 1), with energy centered around that given by Eq. J5J. 

With increasing B concentration, the electron has a 
higher chance to scatter from dimer and trimer impuri- 
ties. For a single dimer impurity, a straightforward but 
cumbersome calculation yields the transmission 
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where N 2 = 
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t 2 AB )(W - 2t AA coskf + t 
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^aa^-bb + ^ab^ 2 ~ 2Wt AA t AB cos k. In general, there 
are more energy values satisfying T 2 (E) = 1, but they 
are different than T\(E) = 1, in general. Consequently, 
although short-range off-diagonal correlation leads to ex- 
tended states for finite systems (localization length l c 
larger than the system size L), these states are not ex- 
tended states in the thermodynamic limit, L — > 00. 

Figure 1 also illustrates that for special local correla- 
tions there is an extended state even for infinite systems. 
One can easily verify that when t AB — t<j = \Jt AA t BB 
(we call this the "golden correlation" in off-diagonal pa- 
rameters), the condition Ti(E) = T 2 (E) = 1 can be 
satisfied. This implies that for the peculiar golden cor- 
relation t(j, single and dimer impurities are essentially 
transparent at this energy. The question remains of 
how general is this result, that is, how about trimers or 
more general impurities? Instead of calculating T m (E) 
for m impurities, we prove the existence of an extended 
state by explicit construction. It is not difficult to check 
that under the condition t AB — to, the state a n e ikn , 
with a n = 1, for e n = e A , and a n = ^t AA /t BB , 
for e n = e B , is indeed an extended state with energy 
E = E c = sa + 2t AA cos k = e B + 2t BB cos k. The physi- 
cal picture for this state is then that the electron propa- 
gates on island A or B in the plane wave form, while the 
golden condition ensures perfect transition from island A 
to island B, and vice versa. One can say that this per- 
fect transmission arises from the cancellation of backscat- 
tered waves produced by the subtle tuning of off-diagonal 
correlations. We find a state with unit transmission co- 
efficient as that shown by curve III in Fig. 1(a), even 
for high concentration of impurities, although this "reso- 
nance" becomes sharper for high impurity concentrations 
(see Fig. 1(b)). This is the first example of an extended 
state in the thermodynamic limit in a random ID system 
with short-range off-diagonal correlations (but no corre- 
lation in onsite energies). Notice also that T(E C ) = 1 



even for a system with 50% disorder, as shown in Fig. 
1(b). 

Under the golden correlation condition t AB — to, the 
extended state satisfies 2 cos k = (e B — e A ) / (t AA — t BB ), 
which can be met only when \e B — e A \ < 2\t AA — t BB \, 
resulting in an interesting effect. Usually in the presence 
of only diagonal or off-diagonal disorder, the larger the 
disorder is, the poorer is the transport. The situation 
is quite different for correlated diagonal and off-diagonal 
disorder. To obtain an extended state in the presence of 
the diagonal difference W — e B — e A , the difference be- 
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FIG. 1: (a) Transmission coefficients vs. energy; W = eb — 
sa ~ 1- Curve I (empty symbols) is typical for system with lo- 
cal correlation; hopping constants here tAA = 1, tBB = 1.73, 
tAB = (tAA + tBB)/2 — 1.36. Curve II (solid line near zero) 
is for Anderson limit, with all hopping constants equal (= 1). 
Curve III (solid symbols) is for system with "golden correla- 
tion", tAB = to = VtAAtBB = 1.316, with unit transmission 
at E ~ —1.9. Concentration of eb site is 0.1 in all three 
curves, (b) Same as in (a), but with concentration of eb at 
0.5. Inset in (a): Localization length 1(E) vs. (E — E c ) for 
system with golden correlation. Slope of fitted line is 2. Con- 
centration of eb site is 0.5. Inset in (b): Localization length 
at critical energy l c = l(E c ) vs. D = (\tAB — £g|), where 
to = VtAAtBB is the golden condition. Concentration of e is 
0.5. 
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tween tAA and tsB has to be large enough, i.e., one needs 
the correlated off-diagonal disorder to be large. This is 
contrary to expectations. 

Notice also that for fixed tAA , and ts b , there is a crit- 
ical onsite difference W = 2\tAA — iss| = 2A. From 
the time evolution of a particle initially placed at a ran- 
domly chosen site (not show here), we find that when 
W < 2 A, the mean square displacement in time r is 
(a; 2 ) ~ t 3 / 2 , and then it is in a superdiffusive phase. 
In contrast, when W = 2A, the system is in a diffusive 
phase, (x 2 ) ~ t; and for W > 2A, the mean square dis- 
placement is bounded. This transition is similar to that 
in the random dimer model (RDM), although with dif- 
ferent characteristics. 2 In RDM, the transition occurs at 
W = 2tAA (all t the same). In our case, the condition is 
related to the difference between hopping constants, and 
not the hopping constants themselves. It is interesting 
that for tAA = tsB < W/2, there are extended states in 
RDM, but no extended state in our model. 

We study the localization length 1(E) for states near 
the critical energy E c in the Fig. 1(a) inset. We find that 
1(E) oc (E — E c )~ 2 for states near E c . The number of 
extended states for a system of length L (i.e., 1(E) > L) is 
related to 5k oc E—E c cx L -1 / 2 , where near E c , E = E c + 
ASk. The number of extended states is then 5k/(l/L) = 
L 1 / 2 , a sizable number, just as in the RDMiiS 

The long time behavior of the system is determined 
by a critical exponent. One can show the relation be- 
tween two exponents 9 and 7, defined by (a; 2 ) ~ t 6 , 
and 1(E) ~ \E — E C \~ J . For short times, the elec- 
tron has ballistic behavior, since it has not sampled yet 
the disorder potential, so that (a; 2 ) ~ (vt) 2 . For long 
time, however, (x 2 ) ~ l 2 (E), for an electron with energy 
E. We can then write the mean square displacement 

as (x 2 ) = J dEp(E)(vr) 2 f (j^^j . wher e p(E) is the 

density of states, and we surmise the scaling function 
f(x) — > 1, as x — > 0, and f(x) — > 1/x 2 , asi-> 00. From 
this, one obtains (x 2 ) ~ r 2-1 / 7 , for long times, so that 
= 2-1/7. When 7 = 2, as in Fig. 1(a) inset, 9 = | 
(superdiffusive regime); while when 7 = 1,9 = 1 (dif- 
fusive). There is perfect agreement with our numerical 
calculations. 

It is natural to expect that in many systems there ex- 
ist correlations between diagonal and off-diagonal disor- 
dered parameters. However the golden correlation con- 
dition is not necessarily satisfied, and it is important 
to see how the transport properties change when a sys- 
tem deviates from this. The inset in 1(b) shows that 
l(E c ) oc (tAB — tc)~ 2 , so that to obtain extended states, 
we need |£ab — ta\ < L^ 1 / 2 . As long as this condition 
is met, effective conduction channels are opened by the 
off-diagonal correlations in the disordered system. 

Our predictions could be verified experimentally in 
systems with access to varying degree of disorder and 
structural correlation, such as model semiconductor su- 
perlattices (SLs)^ Consider a SL with quantum wells 
of two different widths cIa and cIb, distributed randomly 



in the structure. The barriers between wells have the 

same height U (given by the material composition) and 

width bA (or bs) if the barrier is between two alike wells 

of width d,A (or ds), and otherwise have width be- An 

estimate of the hopping constant between two quantum 

wells with width and da, separated by a barrier of 

2*2 

width b and height U is t = m g V^l^_r exp(— where 
s = v2mt/ /h. By tuning parameters, the golden condi- 
tion can be attained. Figure 2 shows the transmission 
for different systems calculated from a Kronig-Penncy 
model of the SL. Curve A is for a system satisfying the 
golden condition, as estimated from the expression above, 
while curves B and C are results away from the condi- 
tion. The discussion above for the tight-binding model 
suggests that transport would indeed be better for system 
in curve A, even as the barrier between different quantum 
wells (curve B) is thinner. We emphasize that Fig. 2 is 
obtained from a Kronig-Penney model of the structure, 
so that hopping amplitudes go far beyond nearest neigh- 
bors, and the golden condition is likely much more in- 
volved than in the tight-binding model. The golden con- 
dition for curve A was not optimized, but just estimated 
from the relation above, and the difference between these 
curves is remarkable^ 

As discussed before, our studies have direct application 
to models of transport in DNA in the literature>ii*i^ For a 
typical DNA molecule the base pair sequence may be es- 
sentially random, such as in A-DNA. However, the chem- 
ical structure determines the local correlation between 
onsite energies and hopping constant via the 7r— orbital 
overlap. In order to explore how the local correlation 
changes transport, we compare the I-V curves of dif- 
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FIG. 2: Transmission for SL with 100 randomly distributed 
quantum wells of two types, width 2.6 (type a well) and 2.9nm 
(type b). Barriers between same a (or b) wells have 3.6nm 
(2.4nm) width. Other barrier width is 3.0nm (curve A, ~ 
golden condition), 2.0nm (curve B), and 3.8nm (C). All bar- 
riers have height 0.3 eV. T(E) averaged over 600 different 
disorder configurations. Concentration of b wells is 0.5. 
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ferent systems, obtained using the Landauer-Biittiker 
formalism^ / = (2e/h) f dET(E)[f L (E) - f R (E)], 
where f L /R.(E) = {exp[E - HL/nJk B T] + l} -1 is the 
Fermi function. We choose [Ib = Ep + (1 — K)eV, and 
/ifl = Ep — neV, where Ep is the equilibrium Fermi 
energy, V is the applied voltage, and k is a parameter 
describing the possible asymmetry of contact to leads, 
chosen here as k = l/3i^& We assume that the DNA is 
attached to ideal leads described by a metal with band- 
width 1.2eV. The hopping constant between leads and 
DNA chain is chosen to be ~ t^s/lO — O.OleV, reflect- 
ing a relatively poor contact. We use two different sets of 
parameters: tAA = — 0.0695eV, tss = — 0.f409eV, and 
tAB = (tAA+tBB)/2, curve A in Fig. 3, describe a realis- 
tic molecule, as the values are obtained from microscopic 
calculations; 10 tAA = Ibb = tAB = — 0.f403eV, curve 
B, simulates an uncorrelated system, i.e., the Anderson 
limit. 

We can see from Fig. 3 that the current in the system 
with local correlation (curve A) is overall much larger 
than in the system without correlation, even though the 



hopping constant is larger in B. There is in fact no con- 
ductance over the entire bias range for curve B, with 
no correlation in the hopping constants. The message 
of these results is that even for DNA with random se- 
quences, such as A-DNA, "good" transport is possible 
due to the effective conduction channels opened by struc- 
tural correlations. Notice that tAB in curve A does not 
satisfy the golden condition (~ — 0.099eV) by about 4%, 
and yet, there is significant current amplitude for finite 
biases. In contrast to the conducting states in polymers, 
which arise from the correlation in local energies, the con- 
ducting states here have to do with correlation in hopping 
amplitudes. It is clear that the backbone may change the 
local correlations. We may conclude that changes in local 
correlation will lead to changes in the I-V features, which 
may in fact be an ingredient in recent experiments, espe- 
cially if chemical changes affect the molecule structured 
We acknowledge support from DOE grant no. DE- 
FG02-91ER45334, NSF NIRT grant no. 0103034, and the 
CMSS Program, and discussions with the NIRT group at 
Ohio U. 
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is still oc L 1/2 . 
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FIG. 3: I-V curves for a random base pair sequence (i.e., ran- 
dom onsite energies). Curve A is for model of A-DNA with 
realistic local correlation in the hopping amplitudes. Curve B 
is for random diagonal Anderson model with hopping ampli- 
tudes set equal. Size of systems is 562; temperature in Fermi 
broadening is 300K. 
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Despite averaging over 600 configurations, curve A in Fig. 
2 exhibits large oscillations for E ~ E c , reflecting the sensi- 
tivity of these states to boundaries and disorder, which pro- 



duces poor self-averaging and arises from their extended 
nature. 



